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Abstract: The next-to- leading order (NLO) QCD corrections to the production and decay 
rate of a Higgs boson are computed within the framework of the Minimal Supersymmetric 
Standard Model (MSSM). The calculation is based on an effective theory for light and in- 
termediate mass Higgs bosons. We provide a Fortran routine for the numerical evaluation 
of the coefficient function. For most of the MSSM parameter space, the relative size of the 
NLO corrections is typically of the order of 5% smaller than the Standard Model value. 
We exemplify the numerical results for two scenarios: the benchmark point SPS la, and a 
parameter region where the gluon-Higgs coupling at leading order is very small due to a 
cancellation of the squark and quark contributions. 
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1. Introduction 



The most important Higgs production cross sections at the Tevatron and the Large Hadron 
Collider (LHC) are known quite reliably within the Standard Model framework. Weak 
boson fusion as well as associated production of Higgs bosons with a top quark pair are 
known with NLO QCD accuracy (see Refs. [1,2] and [3-5], respectively), while gluon fusion 
and Higgs-Strahlung are even known through NNLO in QCD (see Refs. [6-15] and Ref. [16]). 
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At this level of accuracy, also electro- weak corrections can be important. They have been 
evaluated to first order for the Higgs-Strahlung process [17] (a combined analysis of QCD 
and electro- weak effects can be found in [18]). For the gluon fusion process, only very 
recently the full set of first order electro- weak effects has been completed (for < 2Myy): 
The contribution from light fermion loops at two-loop order was evaluated in Ref. [19], the 
top quark induced effects were calculated in Ref. [20]. The terms of order Gprnf have been 
known for ten years [21], and additional QCD effects of order a s GFm^ can be extracted 
from the result of Ref. [22] . 

One remarkable fact about the gluon fusion process is that it has no tree-level contribution. 
It is therefore sensitive to new particles that can mediate the gluon-Higgs coupling. In the 
Standard Model, this coupling is dominated by a top quark loop, with a small contribution 
also from bottom quarks. In the MSSM, the bottom quark contribution can be enhanced 
for large values of tan/3. Their NLO effects have been evaluated in Refs. [8,9]. 

In this paper, we consider the gluon-Higgs coupling mediated by quarks and squarks at NLO 
and its effects on the hadronic production and decay of a light or intermediated mass scalar 
Higgs boson. A preliminary study of these effects has been published in Ref. [23], where, 
however, only a very restricted parameter range of the MSSM has been used. In particular, 
mixing in the stop sector was neglected. In this work we dismiss these constraints. In 
combination with Refs. [8, 9] for the bottom loops, we thus provide the NLO result for 
almost all of the SUSY parameter space. The only restrictions are in the region where 
sbottom-effects become important. However, this only happens when the sbottom mixing 
angle is 6b « 45° and both the mass splitting between b\ and 62 as weu as t an P are very 
large. 

The outline of the paper is as follows. In Sect. ||, we construct the effective Lagrangian for 
the gluon-Higgs interaction by integrating out the top quark as well as all supersymmetric 
particles of the MSSM. The Lagrangian at leading order in 1/M (M G {rat, m h m g}) 
contains only one operator. We evaluate its universal Wilson coefficient through NLO. In 
Sect. |3|, we first study the behavior of the NLO Wilson coefficient itself, and subsequently 
use it to evaluate the hadronic decay and production rate of the light, CP-even Higgs boson 
of the MSSM. To this aim, we choose two specific sets of SUSY parameters: The first one is 
the benchmark SPSla, defined in Ref. [24]. The second one is similar to the "gluophobic 
Higgs" scenario of Ref. [25], where the squark and the quark contributions to the gluon- 
Higgs coupling nearly cancel each other [26]. In Sect.Q, we present our conclusions. App.[A| 
and [B] collect the Feynman rules, counter terms, and decoupling constants that we used 
in our calculation. App.[C] describes the Fortran routine that we provide to evaluate the 
first and second order Wilson coefficient. 
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2. Effective Lagrangian through next-to-leading order 



2.1 Effective Lagrangian 

The effective Lagrangian is constructed from the full MSSM Lagrangian by integrating out 
all SUSY partners and the top quark. The field content of the effective theory is thus the 
same as when starting from the Standard Model Lagrangian. Therefore, also the effective 
Lagrangian has the same form. It is given by 

Adf = ~ cf of, of = \g1^g^ , (2.1) 

where G® „„ is the bare gluonic field strength tensor, v ~ 246 GeV, and Cf is the matching 
coefficient to the full theory. For the sake of simplicity of the discussion, we focus on 
the light neutral Higgs, denoted h, in this paper. The translation of the formulas to the 
heavy neutral Higgs is straightforward, but the validity of the effective theory approach of 
Eq. ([O]) has to be carefully checked in this case. 

For the Standard Model, the two- loop a 2 s corrections for C\{a s ) have been calculated in 
Ref. [6,7], the a 3 and terms in Ref. [27,28] and Ref. [29], respectively. Furthermore, 
as mentioned in the introduction, electroweak corrections of order Gprn 2 and a s GFm 2 
have been evaluated Refs. [21] and [22], respectively. For the MSSM, the two-loop QCD 
corrections are known in the case of zero squark mixing [23]. 

The QCD renormalization of Cf and Of is discussed in Refs. [29,30] and is given by 

C 1 = Z{ 1 1 Cf, 1 = Z u Of, with Z u =(l-—-) , (2-2) 

where 

P(a s ) = -(^f Po + O(a 3 s ) (2.3) 

is the /3-function of standard {ni = 5)-flavor QCD, with (3q = 11/4 — n;/6. Note that here 
and in what follows, a s denotes the strong coupling constant in standard five- flavor QCD; 
it is a function of the renormalization scale //r: 

a 9 = a®(jj, R ). (2.4) 



In this paper we calculate C\ in the MSSM through a 2 , i.e., we will evaluate the coefficients 



cf 1 and Cj defined as 



cS 0) + + 0{a\ 

7T 



(2.5) 



This will allow us to compute the NLO approximation to the hadronic production and 
decay rate of a CP-even Higgs boson. Following the argumentation of Ref. [31], we can 
even derive a fairly accurate estimate of the NNLO production cross section in this model. 
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Several methods to compute the coefficient function C\ are described in Ref. [32]. Here we 
follow the most direct one which is based on the relation 



P1=P2=0 



(2.6) 



T® u (pi,p2) is the 1-particle-irreducible vertex function of two gluons in a color-singlet state 
(incoming momenta p\, P2) and a Higgs boson in the full theory. is the decoupling 
constant that relates the gluon field in the full and the effective theory (details can be 
found in Ref. [32]). It can be computed from the gluon propagator in the full theory U^(p) 
through 



S3 



i + nf(p = o). 



The result for ^ is given in Eq. ( B.lCj ). 
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Figure 1: Diagrams contributing to the effective ggh coupling in the MSSM. 

Sample diagrams corresponding to T® u are shown in Fig. [l]. We may distinguish three 
different types: 



1. pure top contributions, e.g. Fig. [j](a) and (d) 

2. pure stop contributions, e.g. Fig. [j](b), (c), and (e) 

3. mixed top/stop/gluino contributions, e.g. Fig.[l](f)-(i) 
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The pure top quark contributions are separately finite (they correspond to the Standard 
Model terms), while the pure stop and the mixed contributions each develop ultra-violet 
poles that cancel in their sum (after taking into account the proper counter terms). Appli- 
cation of Eq. ( |2.6| ) leads to one- and two-loop integrals with vanishing external momenta. 
They can be evaluated in closed form using the algorithm of Davydychev and Tausk [33] . 



Details will be given in Sect. 2.3 



2.2 Leading order coefficient function 

The LO approximation of the coefficient function is obtained from the one-loop diagrams 
of Fig. 0(a), (b), and (c). The result is 



cos a 




cos 29 t 



2m? 
ti 



m~ 



mi 



m~ 



ta 



(2.8) 



Here and in the following we assume all masses in the on-shell scheme, if not stated oth- 
erwise (mt, mf j and m g are the top, stop (z = 1,2) and gluino mass, respectively), a is 
the mixing angle between the weak and the mass eigenstates of the neutral scalar Higgs 
bosons, tan (3 is the ratio of the vacuum expectation values of the two Higgs doublets, 
Msusy is the coefficient of the bilinear Higgs term in the MSSM superpotential, cf w and 
cf w are defined in Eq. ( |A.7p , and the mixing angle 9t between the helicity (tL,tu) and 



mass eigenstates (ii, £2) of the top squarks is defined in Sect. A..1 . For more details on the 
MSSM parameters, see Ref. [34], for example, q] and ej°~ are the top- and stop-loop con- 
tributions, respectively. Note that the latter are not necessarily suppressed for large stop 
masses. In fact, due to the term cx sin 2 2#t, they can be dominant for large stop mixing 



and large mass splitting between t\ and t<i [26]. Using Eq. ( A. 2 ), this term can be written 
as 



sin 2 20, 



2 

m~ 
ti 

2m? 

ta 



m- 

ta 

2m| 



mi 



4m? m? 



(2.9) 



il ta 



implying that the squark effects are large when X t = A t — ^susy cot is large. In addition, 
Eq. ( ^q ) shows that large effects can also arise from large values of /isusY- 

The question arises if cancellations between the quark and squark contributions occur also 
when radiative corrections are included, or if the regions where this occurs are significantly 
different from the LO prediction. We will discuss this issue for a specific example in Sect.|3|. 
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2.3 Next-to- leading order coefficient function 



A major difference between the Standard Model and the SUSY calculation for C\ is the 
occurrence of more than one mass scale in SUSY; this leads to expressions for C\ that are 
much more complicated and unhandy as compared to the Standard Model result. The 
latter depends only on the top quark mass and thus involves only constants and logarithms 
of the form ln(/x^/m|), where is the renormalization scale. In fact, let us recall the 
expression in the Standard Model [27,28]: 



^,sm _ a s 
61 ~~3^ 



11 a s 
1 + — — + 

4 7T 



2777 19 



Mr 



288 16 



+ — In ^ - ni 



mi 



96 3 m? 



7T 



(2.10) 



where for convenience we also displayed the NNLO result. 

The calculation of in the MSSM leads to two-loop integrals with up to three different 
masses mi, m.2, m 3 £ {m t , , , m^} (integrals with four different masses can be trans- 
formed to integrals with three different masses by simple partial fractioning) . Davydychev 
and Tausk have provided an algorithm for their analytic evaluation [33]. It allows one to 
express the integrals through the function 



4>(mi,m 2 ,m 3 ) = (m 3 A) z 




mi + m2 < m 3 , 
mi + m,2 > m 3 , 



(2.11) 



where 

&i(x,y) 



2 In 



-2Li 2 
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-2Li 2 
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In x In y 
1 2 
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(2.12) 



and 



+ C1 2 



Cl 2 



2 arccos 



-1 + x + y 



2 arccos 



1 + x — y 
2y^ 



xy 
+ C1 2 



2 arccos 



1 - x + y 
2y^ 



A = yj(l-x-y) 2 -4xy, 



mi 



x = 



2 

m\ 



m 



2 ' 



rn 



2 • 



Li2(x) is the standard dilogarithm and Cl 2 (x) is Clausen's integral function, 



Li 2 (x) = - f 1 
Jo 



dt 



ln(l - xt) 



CI 



Jo 



'0 t 

Note that <&(mi, m.2, m 3 ) is symmetric in mi,m 2 and m 3 . 



dt ln|2sin(t/2)| . 



(2.13) 



(2.14) 



(2.15) 



To regulate the ultra-violet divergences of the loop integrals, we use Dimensional Reduction 
(DRED). This is realized by evaluating all Dirac traces and Lorentz contractions in four, 
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and all loop integrals in d = (4 — 2e) space-time dimensions [35,36]. The external projector 
denned in Eq. ( |2.6] ) is taken in d dimensions. Renormalization is done as explained in 
App-i- 

Note that even though we only keep the top- and stop-Higgs couplings different from zero, 
squarks of other flavors (b, c etc.) may enter the NLO calculation through the four-squark 
vertex listed in App. A.2 . Typical diagrams are shown in Fig. [2]. However, diagrams like 



Fig.^(a) vanish due to their color factor, and the diagrams like the one in Fig.|2|(b) add 
up to zero. 



q 



(a) (b) 
Figure 2: Diagrams involving q, with q S {u, d, c, s, t, b}. 

The full result for c± in the MSSM is too long to be quoted here. 1 Instead, we provide 
a Fortran code, named evalcsusy.f [38], that allows for a numerical evaluation of the 
coefficient function and can be combined with other programs quite easily using an SLHA- 
like interface (SLHA: SUSY Les Houches Accord [39]). For details, see App.|C|. 

As a check of our result, we also calculated the diagrams by means of asymptotic expan- 
sions, using the program EXP [40]. It allows us to evaluate an approximate result for C±, 
provided that there is a certain hierarchy among the masses mt, m^, , and m§. The 
approximation, however, will only be valid within the radius of convergence of the specific 
series, so that we will not make use of it in our phenomenological analyses below. Needless 
to say that the expansion of the analytic result expressed through Eq. ( [2.11 ) agrees with 
the corresponding result obtained through asymptotic expansions. 

As another check, we reproduced the results of an earlier publication of ours [23] which 
was obtained by asymptotic expansions and in a very simplifying limit. 



3. Results 

3.1 Results for d 

In order to get an impression about the typical size of the corrections we consider two 
scenarios. First, we look at the behavior of cf^ and at and along a "Snowmass Point 
and Slope" (SPS) [24]. In the second case we consider a particular region of the parameter 
space where C\ shows large deviations from its Standard Model value. 

X A11 algebraic manipulations were done with the help of FORM [37]. 
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To be specific, let us assume an mSUGRA scenario 2 with the five parameters mo, m.1/2, Aq, 
tan/3, and sign(/isusY)- In addition, we define the following Standard Model parameters: 



M z = 91.1876 GeV, m 6 = 4.2GeV, m t = 178 GeV 
«qed( M ^) = 127.934 , G F = 1.16637 • 10 -5 GeV" 2 



m T = 1.777 GeV, 
a s (M z ) =0.118, 



(3.1) 



where Mz is the Z boson mass, mj, = mb(fhb) is the scale-invariant MS value of the bottom 
quark mass, mt is the pole mass of the top quark and m T is the mass of the r lepton. 
&qed(Mz) is the running electromagnetic coupling at Mz, Gf is the Fermi constant and 
a s the strong coupling. As discussed in App.[B|, we further need to define the (arbitrary) 
scale qo (see Eq. flB.7|) ) that enters the renormalization constant of the stop mixing angle 
0t, as well as the usual renormalization scale /^r,. As our default values we adopt 



- (m k + m h ) 



(3.2) 



Mh is the mass of the light CP-even Higgs boson. Both of these choices are generally 
considered to be typical values which avoid the explicit occurrence of large logarithmic 
corrections. 

As already mentioned above and explained in detail in App. |C], the input and output 
files of the program evalcsusy.f follow the SLHA conventions [39]. Among the various 
SUSY-spectrum calculators which are currently available [41-44] (see, e.g., Ref. [45] for a 
comparison), only SoftSusy [42] and SPheno [44] support the SLHA conventions for both 
in- and output. Therefore, we will use these two generators in our analysis. For our 
applications they provide almost identical results. 





-0.002 
-0.004 
-0.006 
-0.008 
-0.01 
-0.012 



-0.014 £ 
-0.016 



SUSY, NLO 
SUSY, LO 
SM, NLO 
SM, LO 



100 150 200 250 300 350 400 450 500 




300 400 500 600 700 800 900 



m 1/2 [GeV] 



(a) 



(b) 



m r [GeV] 

l 2 



Figure 3: LO (dashed) and NLO (solid) result for C\\ (a) as a function of for the SPS la 

scenario; (b) as a function of m t - for the scenario defined in Eq. (3.8). The other parameters are 



fixed as described in the main text. The corresponding Standard Model results are shown as thin 
lines. 



2 Lypical GMSB and AMSB scenarios as defined by SPS 7 and SPS 9 [24] give qualitatively similar results. 



-8- 



SPS la is defined through the following input parameters: 



m = 100 GeV , 
m x /2 = 250 GeV , 

A Q = -100 GeV, 
tan/3 = 10, 

sign(^susY) = +1 • (3.3) 

Using Sof tSusy or SPheno to derive the low energy parameters that enter our result and 
passing them to evalcsusy.f as input, one finds 

cS 0) «1.03, c[ 1} «2.44, (3.4) 
which is close to the Standard Model values 



c 



S 0) = l, c « = ±± = 2.75. (3.5) 



11 



The slope corresponding to SPS la is given by 

rriQ = —A = 0.4m^ 2 , ^1/2 varies . (3-6) 

The dependence of C\ along this slope is shown as thick lines in Fig. || (a) at LO (dashed) 
and NLO (solid). One observes a moderate increase in magnitude of about 8% when going 
from LO to NLO. The thin lines correspond to the Standard Model results. The small 
variation of the latter is due to their dependence on Mh through a s (Mh). For completeness, 
let us remark that the masses that enter our calculation change monotonously within the 
following ranges when going from mi/2 = 100 GeV to mi/2 = 500 GeV: 



101 GeV <M h < 118 GeV, 
176 GeV < m h < 784 GeV, 
330 GeV < m h < 1019 GeV , 
268 GeV < m g < 1158 GeV. 



(3.7) 



The dependence of m^ , tut, and m g on mi/ 2 is almost linear. The dependence of on 
rrii/2 is shown in Fig.||(a). 

In a second example, we consider a case where the LO squark and quark contributions to 
the gluon-Higgs coupling largely cancel each other [26]. Thus, we do not refer to any SUSY 
breaking scenario, but directly choose the following low energy parameters: 3 



m t = 178 GeV , m { = 200 GeV , m g = 1000 GeV , M w = 80.425 GeV , 

7T (3.8) 

- . 250 HeV <m,r < QDD HeV . 
4 



tan/3 = 10, a = 0, t = - , 250 GeV < m t ~ 2 < 900 GeV . 



A qualitatively similar benchmark point ( "gluophobic Higgs") has been suggested in Ref. [46]. 
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m 1/2 [GeV] m r [GeV] 

(a) (b) 

Figure 4: Mh as a function of (a) TO1/2 along the slope of SPS la, and (b) for the scenario 
defined in Eq. 



The light Higgs boson mass is determined by the approximate two-loop formula [47] 

M 2 = M 2 +M 2 +M 2 , (3.9) 



with 



2 3GW2 f /m?\ , X 2 / 1 X 2 



h ' aa ° vr 2 n H \ \M 2 J \ M 2 ) \M 2 ) M s \ 4M| 



(3.10) 



In this approximation and for 0t = vr/4, the parameters Xt and M5 are related to the stop 
masses through 

Xt = i H - <) ' M * = \ K + m • (3 - n) 



The variation of within the parameter range of Eq. (3.8) is shown in Fig. 0(b). 



The choice of 9t is motivated by the explicit result for cj° in Eq. (|2.8[) (see also Eq. (|2.9j) ), 
where the prefactor of the last term in the first line becomes maximal for 0t = vr/4. The 
expression in brackets vanishes for = . However, in the limit -C ?7i^ 2 a term 
enhanced by m~ /m~ survives which can dominate the result for cl, . This is shown in 
Fig. || (b) where C\ and Cf M are plotted as a function of m^ 2 . One observes a rather strong 
variation of the one-loop result of C% . It even changes sign close to 900 GeV, where 

the SUSY and the Standard Model contributions cancel each other. A similar behavior is 
observed at NLO, where C\ vanishes for m^ 2 ~ 850 GeV. 
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3.2 Hadronic decay rate 



For our numerical analysis we neglect all bottom and sbottom effects. In particular, the 
direct coupling of a Higgs boson to bottom quarks is not contained in our formulae. In this 
approximation the LO result for the hadronic decay of a light Higgs boson is determined 
through the h —* gg amplitude shown in Fig.|l|(a)-(c). At higher orders, also multi-particle 
final states contribute, such as ggg, gqq, ggqq, etc. (q ^ t). We write 



T(h — * hadrons) = T y q 



(1 + 5 PS ) = F c 



.( 



V3vr 



(1 + 5) 



(3.12) 



with 



M?y/2G F 



8tt 



9tMn) + 



mi 



mi 



2mJ 9 ^ A ^ )+ 2mJ 



4m? 
Ml ' 



1,2. 



(3.13) 



The second equality in Eq. (|3.12| ) illustrates our approach: the exact leading order result 
proportional to Fq is factored out, and the corrections are treated in the effective-theory 
approach of Eq. (^4]) . The quantity <5ps contains the real and virtual corrections associated 
with the operator 0\. The third equality in Eq. (3.12) is obtained by expanding the ratio 
Ci/c^ in terms of a s . The coupling constants g\ and g^^ in Eq. ( |3 . 1 3| ) are defined in 
Eqs. (ggHgjg), and 



where 



A{t) = -t [1 + (1 - t) f(t)} , A(t) = -~t [l-t f(t)} 



(3.14) 



f(t) 



arcsm 



In 



Vt 



VK 



t>l, 
t<l. 



For completeness, we remark that the limits for t — > oo are given by 



lim Ait) = 1 , lim A{t) = - , 



(3.15) 



(3.16) 



and thus, 



M%V2G F 
8tt 



J0),2 
-i I 



for rat-, mx, m„ 3> , 



(3.17) 



where is given in Eq. 



The quantity 5 is expanded in terms of a s as follows: 

2 



+ 



7T 



7T 



6^ + Oi 



(3.18) 
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and similarly for <5pg. The relation between <5pg and 5 is given by 

^>=4'J+4), *»-®+4f® + 4f + (!|Y, (3.19) 



c 



l °i \°i 



with [6,7,27] 



^)- 73 

(2) _ 37631 _ 363 _ 495 2817 & 363 2 ^ 

dps " 96 8 8 Q[6)+ 16 M 2 + 16 ln M% 

( 7189 11,. 5 , . 263, «?, 11 , 2 ujt " 

+ ^ + — C(2 + -C 3) In n 2 

\ 144 2 SU 4 SU 12 M 2 4 M 2 



(3.20) 



/127 1 w , 7 ^ ill 1 , 2 al 

C(2 + — In + — In 2 ^ 

V108 6 v ; 12 M 2 12 M 2 



where n/ = 5 in our case. 

c-[ is not known in the MSSM, thus only the NLO result for Tg can be calculated consistently 
up to now. However, along the lines of Ref. [31], one can argue that the numerical influence 

(2) (2) (2) SM 

of c\ is small at NNLO, and that it is justified to assume c\ = c\ as long as this 
coefficient has not been computed in the MSSM. The motivation behind this procedure is 
two-fold: On the one hand, one reduces the dependence of the final result on the unphysical 
scales (this is more important for the production rate to be discussed below). On the other 
hand, the relative numerical influence of the coefficient through Eq. (|3.19| ) is more 
important at NNLO than at NLO. Therefore, in the following we will give the numerical 

(2) (2) SM 

result for both NLO and the estimated NNLO (by setting c\ = c\ ). In order to 
indicate that this is not the full NNLO result, we denote it by NNLO'. 

Let us in the following discuss the numerical impact of c$ to the hadronic Higgs decay 



rate in the two scenarios discussed in Sect. 3T. In Fig. ||(a) the decay rate is shown as a 



function of m\j2 where the thick dotted, dashed, and solid line correspond to the LO, NLO, 
and NNLO' prediction for the SPS la scenario. For comparison, we show as thin lines the 
corresponding Standard Model results. 

Similarly, Fig. ||(b) shows the decay rate as a function of m^ 2 for the scenario of Eq. d3.§|). 
As expected, around 850 GeV Tg is close to zero. Furthermore, one observes a screening of 
approximately 50% and more for m^ 2 > 600 GeV. Note that according to Fig. ||(b), it is at 
r 600 GeV where Mh assumes its maximal value; the region above rrif ~ 800 GeV, on 



mr 



the other hand, is experimentally excluded due to the lack of a light Higgs signal at LEP. 
3.3 Hadronic production rate 

In the Standard Model, the total cross section as derived from the effective Lagrangian 
of Eq. (|2.1| ) was shown to approximate the full result to better than 3% for Mjj < 2mt if 



analogous to Eq. ( 3.12 ) — the full top mass dependence at LO is factored out [9,48] 
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Figure 5: Hadronic decay rate for the light CP-even Higgs boson h (a) along SPS la, and (b) for 
the scenario of Eq. The thick lines correspond to the SUSY case, the thin lines show the 

Standard Model results for comparison. Solid: NNLO'; dashed: NLO; dotted: LO. 

(see also Ref. [49]). In this case, the hadronic cross section a = a(pp — > h + X) for Higgs 
production can be written as 

2 



(3.21) 



'a. 



Po 



:',ir) 



A (°) {z) + ^± A m {z) + (?±y a^(z) + ... 

7T V 7T / 



where 
eW(z) 



/ dx\ / dx 



ip%(x\) tpj(x 2 ) S 



(«) 



*j'e{g,q,3} 



z = 



M 2 h 



(3.22) 



(^i(a;) is the density of parton i inside the proton, is the Higgs boson mass, s is the 
hadronic center-of-mass (cm.) energy, and 



7T 



Po 



8M| 



(3.23) 



with from Eq. ( |3.13| ). In order to evaluate the LO, NLO, or the NNLO cross section, the 
second line in Eq. ( |3.21| ) has to be truncated after the term A^°\ or A^ 2 \ respectively. 
Furthermore, the parton density functions (PDFs) <pi in Eq. ( p. 22 ) have to be used at the 
appropriate order. 5 This results in different values for ~E^ n \z) and A^ n \z), depending on 
the order that is being considered. The same is true for a s which has to be set in accordance 
with the PDF set. Specifically, we adopt the PDF parameterizations of Ref. [52,53] where 
a s is given by 0.1300, 0.1165, and 0.1153 at LO, NLO, and NNLO, respectively. 

4 Note that in Eq. (31) of Ref. [23] a factor 1/M| is missing. 

5 Only approximate NNLO parton densities are currently available; with the full NNLO splitting functions 
being known analytically now [50,51], this shortcoming is expected to be eliminated in the near future. 
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Note that at this point we choose a different value for a s (Mz) as the one defined in 
Eq. (|3,1|). The latter enters the evaluation of the low energy parameters through SoftSusy 
or SPheno. This may be viewed as an inconsistency, but we find it more natural to have 
the same set of SUSY parameters at the various orders of the calculation. Besides that, the 
spectrum calculators — to our knowledge — do not provide control over the order of the 
evolution equations, and the numerical effects of the value for a s used in Eq. Q3.1| ) on the 
Higgs production cross section are small. 



The LO partonic result is 



S ig 5 jg 5(1 



(3.24) 



The NLO quantity S^j (x) can be derived from the Standard Model expression of Refs. [6,7]: 



6C2 5(1 - x) + 12 

6(x 2 + 1 -x) 2 
1 — x 



ln(l — x) 



1 — x 



12x(-x + x 2 + 2) ln(l - x) 



ln(x) 



11 



(l-xf 



19 
J qq 1 



| (1 + (1 _ l)2)ln _j_ 



(3.25) 



1 + 2x x l , 

3 ' 



x) 



32 
27 



(1 - xf 



(2) 

The expression for E^- is too long to be quoted here. It can be extracted from Refs. [13-15]. 
In analogy to the discussion below Eq. ( 3.20| ), we define an approximate NNLO result by 

(2) (2) SM 
setting c\ = c\ and denote it by NNLO'. 

For convenience, we explicitly give the relation between the A^ n ^ and the E^ 71 ', n = 0, 1, 2: 



A (o) 



<0) 



A ( 1 ) = s( 1 ) + 2^ y s(°) 



A( 2 ) = S( 2 ) + 2^ EW + 



' c (i)y c (2) 

Jo)) +2 JV 



1 



-1 



(3.26) 



E (0) 



The quantities E^ ra ^ are independent of the specific model under consideration. A publicly 
available numerical program for their evaluation is in preparation [54]. 

For illustration of the numerical effects on the total Higgs production cross section in 



gluon fusion, we consider again the two exemplary cases of Sect. 3.1. The thick lines of 
Fig. ^(a) show the LO (dotted), NLO (dashed), and NNLO' (solid) cross section for the 
SPS la scenario, compared to the Standard Model prediction (thin lines). At m 1 / 2 = 
100 GeV, the MSSM values are about 34% (32%) larger than the Standard Model values at 
NLO (NNLO'), but the difference decreases quite rapidly when mj/j increases. However, 
this difference is mostly a LO effect, as can be seen from Fig. || (b) which shows the NLO and 
the NNLO' K- factor in the MSSM and the Standard Model. The difference of the relative 
corrections in both cases is less than 5%. 
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Fig. 0(a), on the other hand, corresponds to the scenario defined in Eq. (3.8). As seen in 
Fig.||(b), C\ vanishes for a certain value of ?7ir , and so does the NLO and the NNLO' 
cross section. Note, however, that this particular value is experimentally excluded because 
it corresponds to a too low Higgs mass (see Fig. |](b)). Nevertheless, for m ( - 2 ~ 600 GeV, 
where is maximal, the cross section is still significantly suppressed with respect to the 
Standard Model. As the K-factor in SUSY tends to be a little smaller than in the Standard 
Model, this suppression becomes even stronger when QCD corrections are included. For 
example, at mj ~ 600 GeV (or alternatively, \X t \ 
0.58 at LO, 0.52 at NLO, and 0.48 at NNLO'. 



900 GeV), the ratio a MSSM /a SM 



is 
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Figure 6: (a) Total cross section at LO (dotted), NLO (dashed), and NNLO' (solid) along the slope 
of SPS la (thick lines). For comparison, also the Standard Model result is shown (thin lines). It 
depends on mi/ 2 due to the variation of M/j with this parameter, (b) NLO and NNLO' K-factor for 
the SPS la scenario (thick lines) and for the Standard Model (thin lines). 
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Figure 7: (a) Cross section and (b) K-factor for the scenario defined in Eq. (L8) (thick lines) and 
for the Standard Model (thin lines). Dotted lines correspond to LO, dashed lines to NLO, and solid 
lines to NNLO'. 
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3.4 Discussion 

In a model where the gluon-Higgs coupling is mediated predominantly by heavy particles, 
it had already been observed that the radiative corrections to the hadronic production and 
decay processes are not very sensitive to the specifics of this coupling [31,55]. This is due 
to the fact that the radiative corrections are dominated by soft gluon effects which do not 
resolve the gluon-Higgs vertex [11-13, 56]. 6 

Aside from this, for typical MSSM benchmark points, even the Wilson coefficient of the 
effective gluon-Higgs coupling itself is numerically rather close to its Standard Model value, 
both at LO and NLO. Only if at least one of the scalar top quarks is relatively light 
{ m t ~ 400 GeV), a significant deviation from the Standard Model result is observed. This 
is because the stop Yukawa coupling is proportional to m\ rather than m~. In combination 
with the loop amplitude of Fig.|](b),(c), this leads to a suppression factor mf/m|. In 
contrast to this, for the top quark contribution there is a cancellation between the Yukawa 
coupling ~ mt and a factor 1/mt from the loop amplitude. 

We pointed out that this suppression of the squark contribution may be compensated by a 
large absolute value of the parameter X t = A t — /^susy cot (3. According to Eq. (|2.9|) , this 
corresponds to a stop mixing angle of the order of 0t = vr/4, and a large mass splitting 
between mi and m^ 2 . However, the value of Xt is crucial for the exact value of the light 
Higgs boson mass M^. This restricts \Xt\ to less than about 3TeV. In Eq. ( |3.8[ ), we have 
chosen a set of low energy parameters which fulfills this condition, but where the Wilson 
coefficient C\ for the gluon-Higgs interaction is very different from its Standard Model 
value (see Fig. || (b)) and leads to a strongly reduced production and decay rate. Since also 
here the QCD corrections tend to be smaller than in the Standard Model, this cancellation 
effect of top and stop contributions is even stronger when QCD corrections are included. 

4. Conclusions 

We have analytically calculated the NLO QCD contribution to the effective gluon-Higgs 
coupling in the MSSM due to the scalar partners of the top quark. Scalar bottom effects are 
generally suppressed by m|/m~ or m^susY /mf and have been neglected. The calculation 
involves Feynman diagrams with three massive particles (gluino, top quark, stop quark) 
which leads to very long analytic expressions for the final result. Therefore, we make it 
available in the form of a Fortran routine, described in App.^. 

The results for the effective coupling were used to evaluate the hadronic Higgs decay rate 
and the production cross section through NLO in QCD, and to derive a NNLO estimate 
of these quantities. The QCD corrections in the MSSM tend to be a bit smaller than 
in the Standard Model. However, this effect is in general below 5%. In regions of the 

6 It is remarkable, however, that the resummation of these soft-gluon effects gives only an effect of about 
6% [56] with respect to the fixed order NNLO result. This is usually interpreted as a sign of the stability of 
the NNLO prediction. 
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MSSM parameter space where the Higgs coupling to gluons is particularly small due to a 
cancellation between the quark and the squark contribution, the reduced K-factor amplifies 
this effect. Nevertheless, even here the K-factor in the Standard Model provides a fairly- 
accurate approximation to the MSSM value. 

We conclude by noting that the methods of our calculation should be immediately appli- 
cable to the photonic production and decay rate of an MSSM Higgs boson, as well as to 
pseudo-scalar Higgs production. Inclusion of sbottom effects is also possible, but requires 
a careful treatment of the bottom threshold in the Feynman diagrams. 

Acknowledgments. The authors are grateful to M. Faisst, S. Heinemeyer, and S. Kraml 
for useful discussions. RVH is supported by Deutsche Forschungsgemeinschaft, contract 
HA 2990/2-1. 
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A. Feynman rules 



In this appendix, we collect the Feynman rules that have been used in our calculation of 
the NLO coefficient function C\. The notation follows closely Ref. [57]. 



A.l Definitions 

In the following, p, k, and p n (n = 1,2,3) denote incoming four-momenta; the various 
indices have the following meaning: 



r, s, t, u 


: color triplet indices 


a, b, c 


: color octet indices 


H,v,p 


: Lorentz indices 


i, j, k, I 


: squark mass eigenstate indices 


A, B 


: flavor indices 



Furthermore, we introduce 



[T a , T b ] = if abc T c , {T a , T b } = —5 ab + d abc T c , 
Pl = ^, Pn 1 + 75 



TV 
S q ■- 



2 ' n ~ 2 ' 
'K q u K q 12 \ _ ( cos0 9 sin0 ? \ (A.l) 
K\ x TZ q 22 ) ~ y-sme q cos6 q ) ' 

S ll S 12 \ _ I COS 20 Q ~ Sm 26 q \ 
S 2l S 22 ) ~ V - Sin 2 °q ~ COS 2 °q ) ' 



9 q G [0, |) is the squark mixing angle defined through 



2m„a„ rn* — Tr?= 
sin 29 q = 2 g % , cos 29 q = -f , (A.2) 



where, by definition, we assume < mg 2 , and 



v a j cot /3 , for <? G {u, c,t} 

X q = A q -» SVSY . for q€{di8jb} • (A.3) 



Aq and //susy are soft SUSY breaking parameters (see Ref. [34], for example). 
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A. 2 Feynman rules used in this calculation 
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ig a T? a y/2(nlP L -n^P R ) 



1i 




-ids {^abSrs + d abc T£ s ) 5, 



rpa rpa cA cB i rpa rpa cA cA s , 
J rs J tii D ij D fci + 1 ru- L ts' :> il' :> kj AB 



(S A = Si 



We only give the values for the couplings to the light Higgs here; the couplings to the 
heavy Higgs can be obtained from the latter through the replacement a — > a + tt/2. The 
top-Higgs coupling is 



9t 



cos a 
sin/3 ' 



(A.4) 



and the stop-Higgs couplings are 



6 i,^ s, a 



(A.5) 



- 20 - 



with 



where 



h,EW 

9t'n 

h,EW 
9t,22 

h,EW 
9tX2 



h.u 

9 t ii 



h.u 
9 t 12 



h,a 

9 t ,u 



h,a 
9t,22 



h,a 
9 t ,12 



cf W cos 2 ft + cf W sin 2 ft, 
cf w sin 2 ft + cf w cos 2 ft, 

h,EW _ 1 / EW „EW\ ■ O o 
9t,21 - ol c 2 ~ c l ) SHI 2ft 



ft) 



h,n _ /^susy cos (a 
~ 922 ~ mt sin 2 /3 
h,n _ A^susy cos(a - (3) 



sin 2ft , 



9 t ,21 

cos a 
sin /3 

cos a 
sin/3 



in 



sin 2 /3 

,2 



cos 2ft , 



TO; 



2 + \ „ * 2 sin 2 2ft 



2m 2 



2 

51 



m 



2m 2 t 



^ sin 2 2ft 



h)Q _ cos a 
5 *' 21 " sin/3 



ra- 



ti 



2m 2 t 



sin 2ft cos 2ft , 



(A.6) 



EW 
c l 



c 2 



sin 6w 



Mi 



HIT 



1 — - sin 2 #iy I sin(a + /3) 



M l 4 • 2 „ • / 

f- - sin 6\y sm(a + p) , 

m+ 3 



(A.7) 



and 



2A% 



V\/2Gf 



u 2 -f uf « 246 GeV, 



(Ai 



with vi, V2 the vacuum expectation values of the two Higgs doublets. In Eq. ( |A.6| ) we 
have already expressed the trilinear couplings of the soft SUSY breaking terms through 
independent parameters: 



A, 



1 2 

m~ — m~ 

t\ 12 



2m t 



sin 2ft + hsvsy cot (5 . 



(A.9) 



The electroweak radiative corrections to this formula need not be considered here. 



B. Renormalization and decoupling constants 



In order to arrive at a finite NLO result, the parameters appearing in the LO coefficient 
function given in Eq. ( |2.g| ) have to be renormalized. This includes the strong coupling 
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constant a s , the top quark mass mt, the top squark masses m^., and the mixing angle 9t, 
whereas the angles a, [5 and 6w are not renormalized, because we consider QCD corrections 
only. 

First, the top quark and the SUSY partners are decoupled from the bare coupling constant 
through the relation 



« B = (C B ) 2 « B , 



(B.l) 



where a® and denote the bare couplings in the full theory and in five-flavor QCD, 
respectively. In the DR scheme, we find 



(C B ) 2 



7T 



\ CA + \ T 



+ L(e))+0(a 



(B.2) 



where Ca = 3, T = 1/2 and 



L(e) = - [2C A L g + T(L h + L~ t2 + 4L t )] 



e 

+ 12 



C A L~ + -T (L| + L? 2 ) + 2T L? + (CU + 3T) Q 2 



In 



Mr 

2 ' 

m~ 



In^ 

m~ n 



(B.3) 



is then renormalized through 

„B 



at 



a 



DR 



a s 1 
1 + — - 

7r e 



12 



-Tni 
3 



(B.4) 



where a]? R denotes the DR expression for the strong coupling constant in QCD with m = 5 
active flavors. We will comment on the transformation to the more familiar MS scheme 
below. 

For the quark and squark masses, we adopt the on-shell scheme, where they are defined 
as the real part of the pole of the corresponding propagator. Furthermore, we define the 
renormalized squark mixing angle by requiring that the non-diagonal two-point function 
(ti<2) vanishes at a certain momentum transfer qo; i.e., the two squarks propagate inde- 
pendently from each other at the scale qo. In practice, qo is chosen to be of the order of the 
squark masses. The counter terms can be found in Ref. [58], for example. For convenience, 
we list them explicitely in our notation. 

In DRED, the relation between the bare and the pole top quark mass reads 



m 



a 



_L = i + Cf ^{ — — -L 



7T 



2c 



2 



m 2 



m 9 

Ami 



m~ n 



2 



8m 



%(1 + Lti) 



+ i ( 1 + 3 ~ dr + 2 (- X ) l d sin2 ^ ) B^(mi,m- 9 ,m k ) 



m 



m 



m t 



(B.5) 
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where Cp = 4/3. The only difference between DRED and Dimensional Regularization 
(DREG) comes from the gluon-exchange diagram which changes the constant "5/4" into 
"1" in the case of DREG. The relation for the top squark mass mi is given by 



m 



B 



1 



m;, 7r 8mf e 



'fi 



O O/O O \ o 

Arrtgmt sin 2#t — 4m^ — 4m t + (m^ — ) sin 2#* 



+ 



3 sin 2 26 t 

4 8 

m? sin 2 26*/ 

J2 

8m? 

El 



1 sin 2 2(9 # 



L 



ti 



m 9 

4m? 

El 



4m? 

*i 



(l+^ 2 ) + 



1 2m g mt sin 20* — m? — m 2 



4m? 

El 



B n (rri~ ,mt,rrig) 



(B.6) 



The corresponding relation for the mass is obtained by interchanging the indices "1" 
and "2" and changing the sign of sin 20*. 

Finally, for the mixing angle we have 



Of = 9 t + 56 t , 



66, 



C F 



a, cos 20* f 4m 8™i + ( m | - m \ ) sin 

~4e 



+ 



7r m~ 

El 

sin26»* 



ra- 



ti 



m 



(1 + L i2 ) - m| (1 + L fl ) + mgmX n (9o> m*, m g ) \ . (B.7) 



fa 



For q 2 < (mi — m.2) 2 , B^ n (q 2 , mi, m-z) is given by 



B^(q 2 ,m 1 ,m 2 ) 



In 



m\m2 m\ — m\ ni2 M + M_ M + + M. 



+ 



In — + 

mi 



In- 



M, - M_ 



(B.« 



with M± = -y/ (mi ± m2) 2 — (? 2 . The analytical expressions for the other kinematical re- 
gions can be derived from this expression by proper analytical continuation. Note that the 
counter terms in the DR scheme are obtained by discarding the finite parts at order a s . 



The decoupling constant entering Eq. ( [2.6D is defined in analogy to Eq. (B.l) via the relation 
of the bare gluon field of the full theory, G®, (i.e. including the top quark and the SUSY 
particles) and the effective theory, G®: 



G, 



^3 n > 



where 



S3 



1 + — 

7T 



+ Ofa 2 



(B.9) 



(B.10) 



with L(e) from Eq. 

a° R is transformed from the DR to the MS scheme through a finite shift [59]; however, 
we found that this shift is canceled by a finite shift in the decoupling constant and 
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the operator renormalization Zn, given in Eq. (|2.2|), Our final result is thus expressed 
in terms of the MS coupling a s for standard five-flavor QCD, on-shell quark and squark 
masses, and the squark mixing angle as defined in Eq. ( B.7|) (the gluino mass is unaffected 
by renormalization at the order considered here). 



C. Description of evalcsusy.f 



In this Appendix we give details on the usage of the Fortran program evalcsusy.f [38]. 
The distribution includes the files 



evalcsusy.f common-slha. f f unctions. f readslha.f slhablocks.f 

as well as the Fortran code for the two-loop results of C\ as defined in Eq. ( |2.5[ ); the latter 
is contained in the directory obj/. For compilation one also needs the CERN libraries 
kernlib and mathlib. 

In a first step, object files are generated with 

> ill -c f unctions. f readslha.f slhablocks.f obj/*.f 

The executable is then obtained by saying 

> ill -o evalcsusy evalcsusy.f f unctions. o readslha.o \ 
slhablocks.o obj/*.o -L'cernlib -v pro kernlib , mathlib ' 



evalcsusy is invoked as 



> ./evalcsusy <infile> <outfile> 



where <inf ile> and <outf ile> are the in- and output file, respectively. Both files obey 
the SUSY Les Houches accord (SLHA) [39] which makes it straightforward to interface 
evalcsusy.f with a spectrum calculator. The basic idea of the SLHA is to group the 
parameters into various blocks which have a uniquely defined structure in order to ensure 
universality. For our process we need some parameters (the precise specification can be 
seen in the example presented below) of the blocks SMINPUTS, MASS, ALPHA, HMIX, STOPMIX 
and MINPAR. In addition, we introduce a new block CREIN specific to evalcsusy.f, where 
the ratio /ir/M^ and the parameter go is defined (see Eq. (|3.2| ) and Sect.|l|). The latter is 
actually composed out of the three quantities q$, goi an d qo2 via the relation 

<?o = Qo + Qoim^ + qo2m i2 . (C.l) 

If qoi (qo2) is not defined, its value is set to zero. 

evalcsusy.f copies the contents of the input file to the output file and appends an addi- 
tional block HGGSUSY. Its structure is as follows: 
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Block HGGSUSY 



101 

102 

103 

201 

202 

1001 

1011 

1022 

1012 

1021 



„(0),SM 
-1 

„(1),SM 
-1 

_(2),SM 
-1 

JP) 
"1 

,(1) 



9t 

gin 

9t,22 

a h 

y*,i2 

9t,21 



It contains the results for the one- and two-loop coefficients (c.f. Eq. (|2.5|)) both for the 
Standard Model (c^' SM and c^' SM ) and the MSSM (cf^ and c^). Furthermore, the three- 
loop Standard Model term is given. In addition, it provides numerical values for the top 
and stop Yukawa couplings g\ and g^j (see Eqs. (|A.4j) -( [Alq) ). 

Let us exemplify the use of evalcsusy.f by considering the SPS la point given in Eq. (|3.3|). 
Defining in addition the Standard Model parameters of Eq. (|3.l| ) in the block SMINPUTS, 
the output of SoftSusy [42] looks as follows (after adding the block CREIN): 



Block CREIN 

6 O.dO # qOc 

61 .5d0 # qOl 

62 .5d0 # q02 

7 l.dO # renormalization scale muR/mh 



Block SMINPUTS # Standard Model inputs 

3 1.18000000e-01 # alpha_s (MZ)MSbar 

4 9.11876000e+01 # MZ(pole) 

6 1.78000000e+02 # Mtop(pole) 



Block MINPAR # SUSY breaking input parameters 



3 


1 


00000000e+01 


# 


tanb 


4 


1 


00000000e+00 


# 


sign(mu) 


1 


1 


00000000e+02 


# 


mO 


2 


2 


50000000e+02 


# 


ml2 


5 


-1 


00000000e+02 


# 


AO 



Block MASS # Mass spectrum 

#PDG code mass particle 



24 8.02591534e+01 # MW 

25 1.12153306e+02 # hO 
1000006 3.97398225e+02 # ~t_l 
1000021 6.11147741e+02 # ~g 
2000006 5.86830420e+02 # ~t_2 

Block alpha # Effective Higgs mixing parameter 
-1. 13348399e-01 # alpha 

Block hmix Q= 4 . 661391312e+02 # Higgs mixing parameters 
1 3.65690378e+02 # mu 

Block stopmix # stop mixing matrix 
1 1 5.34006091e-01 # 0_{11} 

1 2 8.45480617e-01 # 0_{12} 

where only those parameters are displayed which are needed in evalcsusy.f. Using this 
file as input for evalcsusy.f, its contents are copied to the output file, and the block 



Block HGGSUSY 










101 


0. 


10000000E+01 


# 


cSM 


1-loop 


102 


0. 


27500000E+01 


# 


cSM 


2-loop 


103 


0. 


35160026E+01 


# 


cSM 


3-loop 


201 


0. 


. 10343231E+01 


# 


cSUSY 


1-loop 


202 


0. 


24377938E+01 


# 


cSUSY 


2-loop 


1001 


0. 


99853850E+00 


# 




gth 


1011 


-0. 


53035841E+00 


# 




gthll 


1022 


0. 


42680033E+01 


# 




gth22 


1012 


0. 


. 11983135E+01 


# 




gthl2 


1021 


0. 


. 11983135E+01 


# 




gth21 



is added. 
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